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During  the  past  few  decades,  wave  energy  has  received  significant  attention  for  harnessing  ocean  energy. 
Industry  has  proposed  many  technologies  and,  based  on  their  working  principle,  these  technologies 
generally  can  be  categorized  into  oscillating  water  columns,  point  absorbers,  overtopping  systems,  and 
bottom-hinged  systems.  In  particular,  many  researchers  have  focused  on  modeling  the  point  absorber, 
which  is  thought  to  be  the  most  cost-efficient  technology  to  extract  wave  energy.  To  model  such  devices, 
several  modeling  methods  have  been  used  such  as  analytical  methods,  boundary  integral  equation  meth¬ 
ods  and  Navier-Stokes  equation  methods.  The  first  two  are  generally  combined  with  the  use  of  empirical 
solution  to  represent  the  viscous  damping  effect,  while  the  last  one  is  directly  included  in  the  solution.  To 
assist  the  development  of  wave  energy  conversion  (WEC)  technologies,  this  paper  extensively  reviews 
the  methods  for  modeling  point  absorbers. 

©  2012  Published  by  Elsevier  Ltd. 
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1.  Introduction 

While  development  of  modern  wave  energy  converter  dates 
back  to  1799  [1],  the  technology  did  not  receive  worldwide  atten¬ 
tion  until  the  1970s  when  an  oil  crisis  occurred  and  Stephen 
Salter  published  a  notable  paper  about  the  technology  in  Nature 
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in  1974  [2].  In  the  early  1980s,  after  a  significant  drop  in  oil  prices, 
technical  setbacks  and  a  general  lack  of  confidence,  progress  slowed 
in  the  development  of  wave  energy  devices  as  a  commercial  source 
of  electrical  power.  In  the  late  1990s,  awareness  of  the  depletion 
of  traditional  energy  resources  and  the  environmental  impacts  of 
the  large  utilization  of  fossil  fuels  significantly  increased,  thereby 
facilitating  the  development  of  green  energy  resources.  The  devel¬ 
opment  of  the  wave  energy  technology  grew  rapidly,  particularly 
in  oceanic  countries  such  as  Ireland,  Denmark,  Portugal,  the  United 
Kingdom,  and  the  United  States.  Several  pre-commercial  devices 
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were  deployed.  For  example,  a  United  States  company,  Ocean 
Power  Technology,  deployed  one  of  their  1 50  kW  wave  energy  con¬ 
version  (WEC)  systems  in  Scotland  in  2011  [3].  An  Irish  company, 
Wavebob,  tested  a  one-quarter  scale  model  in  Galway  Bay,  Ireland, 
in  2006  [4].  In  Denmark,  the  half  scale  600  kW  Wave  Star  energy 
system  was  deployed  at  Hanstholm  in  2009  [5],  and  a  quarter 
and  a  half  size  model  Wave  Dragon  was  tested  at  Nissum  Bred- 
ning  in  2003  [6].  Furthermore,  international  organizations,  such 
as  the  International  Energy  Agency  and  the  International  Elec¬ 
trotechnical  Commission,  are  heavily  involved  in  the  development 
of  wave  energy  devices.  In  2001,  the  International  Energy  Agency 


established  an  Ocean  Energy  System  Implementation  Agreement  to 
facilitate  the  coordination  of  ocean  energy  studies  between  coun¬ 
tries  [7].  In  2007,  the  International  Electrotechnical  Commission 
established  an  Ocean  Energy  Technical  Committee  to  develop  ocean 
energy  standards  [8]. 

1.1.  Device  design 

To  date,  there  are  more  than  one  hundred  prototypes  of  various 
WEC  systems  [7].  They  can  be  sorted  into  a  few  types.  According  to 
their  energy  conversion  principles,  WEC  systems  include  oscillat- 
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Fig.  1.  An  illustration  of  the  floating-pitching  devices:  (top)  multiple  degree  of  freedom  and  (bottom)  single  degree  of  freedom.  Illustration  by  A1  Hicks,  NREL. 
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ing  water  column,  overtopping  system,  pitch  device,  membrane, 
and  point  absorber.  Most  of  these  can  be  both  bottom-mounted 
and  floating.  A  brief  introduction  on  the  working  principle  of 
different  types  of  WEC  system  is  given  in  this  section,  but  the  focus 
of  this  paper  is  on  point  absorbers. 

The  floating-pitching  device,  whether  it  is  composed  of  a  sin¬ 
gle  body  or  a  number  of  connecting  bodies,  has  rotational  freedom 
(Fig.  1 ).  The  device  converts  wave  energy  from  its  pitching  motion. 
The  principal  axis  for  floating-pitching  devices  is  either  perpen¬ 
dicular  (terminator)  or  parallel  (attenuator)  to  the  wave  direction. 
Among  various  developments  since  the  1970s,  Salter  proposed 
one  of  the  most  significant  inventions  in  1974  [2],  which  became 
the  well-known  Salter’s  Duck  (also  referred  to  as  the  Edinburgh 
Duck).  The  Salter’s  Duck  and  PS  Frog  Mk  5  [9]  are  two  examples 
of  the  terminator,  and  McCabe  Wave  Pump  [10]  and  Pelamis  [11] 
are  the  examples  of  the  attenuator.  A  wave  farm  consisting  of 
750  kW  Pelamis  devices  was  tested  in  northern  Portugal  in  2008 
and  another  farm  was  deployed  in  Orkney,  Scotland,  in  2010. 

The  bottom-hinged  pitching  system  includes  a  paddle,  or  a  flap, 
that  is  connected  to  a  hinge  deflector  on  the  seabed;  the  top  of  the 
device  is  generally  above  the  free  surface  (Fig.  2).  It  is  sometimes 
called  a  wave  surge  converter,  because  it  converts  wave  energy 
from  the  horizontal  movement  of  the  water  particles.  Unlike  the 
floating-pitching  device,  one  end  of  the  bottom-hinged  device  is 
fixed.  That  is,  it  shares  a  similar  working  principle  as  the  float¬ 
ing  terminators,  e.g.,  Salter’s  Duck  in  Fig.  1  (bottom).  Examples  of 
the  bottom-hinged  pitching  system  include  the  swing  mace,  also 
proposed  by  Salter  in  the  1990s  [12],  and  the  Aquamarine  Power 
Oyster  [13],  which  was  connected  to  the  electrical  grid  and  tested 
in  Scotland  in  2009  [14]. 

The  oscillating  water  column  includes  a  special  chamber  with 
a  bidirectional  turbine  inside  (Fig.  3).  One  end  of  the  chamber 


has  an  inlet  that  allows  the  incident  wave  to  enter  and  the  other 
end  contains  the  turbine.  The  device  converts  wave  power  by 
utilizing  the  wave  elevation  to  compress  or  decompress  the  air 
in  a  chamber.  The  compressed  air  goes  through  a  bi-directional 
turbine.  The  turbine  is  a  Wells  turbine  conceived  by  Professor  Alan 
Wells  of  Queen’s  University,  in  Belfast,  in  late  1970.  A  wide  variety 
of  oscillating  water  column  devices  include  the  LIMPET  (shoreline 
system)  [15],  and  Sakata  (integrated  into  breakwater)  [16]. 

The  overtopping  device  includes  a  large  structure  that  embraces 
the  incident  wave  and  an  outlet  with  turbines  inside  the  large  struc¬ 
ture  (Fig.  4).  The  device  converts  wave  power  by  utilizing  the  wave 
overtopping  phenomenon  to  let  the  water  fall  through  the  outlet 
of  the  designed  structure.  When  the  water  falls  through  the  out¬ 
let,  it  passes  one  or  more  turbines  similar  to  a  traditional  hydro 
dam;  the  potential  energy  is  converted  into  electric  power.  The 
design  involves  both  kinematic  energy  and  potential  energy  in  the 
conversion  process.  Examples  include  the  Tapchan  (shoreline  sys¬ 
tem)  [17],  the  Wave  Dragon  (offshore  floating  system)  [18],  and 
the  Seawave  Slot-Cone  Generator  (SSG),  which  was  integrated  into 
a  breakwater  [19]. 

The  membrane  device  includes  a  membrane  structure  and  a 
power  conversion  system  that  can  be  a  turbine,  piezoelectric,  or 
other  system  (Fig.  5).  The  device  converts  wave  power  by  utiliz¬ 
ing  the  dynamic  pressure  change  during  wave  propagation  (e.g., 
Lilypad  [20]). 

A  point  absorber  generally  converts  wave  energy  by  utilizing 
the  heave  motion  (Fig.  6).  Note  that  there  are  few  pitching  point 
absorbers,  such  as  PS  Frog  Mk5  [9],  and  they  are  considered  as 
pitching  devices  in  this  paper.  Typically,  a  point  absorber  is  either 
a  single-body  that  generates  energy  by  reacting  against  a  fixed 
seabed  frame,  or  it  is  a  multiple-body  device  that  generates  energy 
from  the  relative  motion  between  the  two  bodies.  There  are  many 
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Fig.  2.  An  illustration  of  the  bottom-hinged  system.  Illustration  by  A1  Hicks,  NREL. 
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Fig.  3.  An  illustration  of  the  oscillating  water  column.  Illustration  by  A1  Hicks,  NREL. 


Fig.  4.  An  illustration  of  the  overtopping  device.  Illustration  by  A1  Hicks,  NREL. 


4356 


Y.  Li,  Y.-H.  Yu  /  Renewable  and  Sustainable  Energy  Reviews  16  (2012)  4352-4364 


Fig.  5.  An  illustration  of  the  membrane  device.  Illustration  by  A1  Hicks,  NREL. 


popular  devices  such  as  the  OPT  PowerBuoy  [16],  Wavebob  [21], 
and  Inter  Project  Service  buoy  [22-24].  Since  the  point  absorber  is 
the  focus  of  the  paper,  its  principles  are  described  in  greater  detail 
in  Section  2. 

2 .2.  Objective  of  this  paper 

Several  literature  reviews  of  WEC  devices  have  been  published 
providing  information  about  various  aspects  of  the  technology.  For 
example,  Evans  reviewed  the  analytical  derivation  of  the  absorber’s 
motion  [25].  Falcao  provide  a  comprehensive  overview  of  the  his- 
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Fig.  6.  An  illustration  of  the  point  absorber.  Illustration  by  A1  Hicks,  NREL. 


tory  and  status  of  the  development  of  WEC  systems  [26].  Sarmento 
discussed  the  non-technical  aspects  of  the  technology  during  its 
commercialization  path  [27].  The  U.S.  Department  of  Energy  (DOE) 
developed  a  database  to  document  various  technology  character¬ 
istics  [28].  Some  concise  reviews  reporting  interim  progress  were 
written  by  several  researchers,  such  as  Falnes  [29],  on  theoretical 
limits  [29],  and  Khan  et  al.  [30],  on  interconnection  issues.  There 
also  are  introductory  textbooks  on  this  subject  such  as  Wave  Energy 
Conversion  by  Brooke  [31  ],  Ocean  Wave  Energy  by  Cruz  [32],  Ocean 
Wave  Energy  Conversion  by  McCormick  [33],  and  Ocean  Waves  and 
Oscillating  Systems  by  Falnes  [34]. 

However,  to  date,  there  has  been  no  systematic  review  of  mod¬ 
eling  methods.  The  rapid  development  has  not  led  the  industry  to  a 
commercial  stage  yet.  From  a  pre-commercial  stage  to  a  full-scale 
commercial  stage,  a  good  understanding  of  device  performance  and 
reliability  is  required.  This  requirement  has  not  been  met  yet.  In 
addition  to  these  non-technology  barriers,  it  is  believed  that  a  lack 
of  understanding  and  quantification  of  the  device  behavior  is  one 
of  the  main  reasons  industry  development  has  been  delayed. 

During  the  development  of  wave  energy  technology,  many 
methods  for  modeling  WECs  were  developed  during  the  technol¬ 
ogy’s  evolution.  To  understand  point  absorber  behavior,  precise 
modeling  tools  are  needed  to  simulate  the  dynamics  of  the  device, 
along  with  a  well-defined  combination  of  modeling  methods.  Tech¬ 
nology  developers,  engineers,  and  researchers  need  a  guideline 
of  which  methods  to  employ  for  a  specific  purpose.  To  facilitate 
the  development  of  the  technology  and  help  to  model  the  device 
behavior  more  cost-effectively,  we  conducted  a  systematic  review 
of  various  methods  for  modeling  WECs.  The  reviews  are  summa¬ 
rized  in  this  paper.  Since  there  are  a  great  number  of  device  types, 
we  will  focus  on  the  point  absorber  as  an  example.  Additionally,  we 
focus  on  the  hydrodynamic  modeling  methods  in  this  paper.  Other 
modeling  methods,  such  as  electrical  or  control  modeling,  are  not 
discussed,  although  they  are  important  as  well. 

Specifically,  this  paper  summarizes  the  numerical  methods  used 
for  simulating  point  absorbers.  After  presenting  the  working  prin¬ 
ciples  of  the  point  absorber,  the  paper  details  the  existing  modeling 
work  that  uses  analytical  approaches,  empirical  methods,  and 
numerical  methods.  At  the  end,  we  compare  the  advantages  and 
disadvantages  of  different  hydrodynamic  modeling  methods. 
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Fig.  7.  The  translation  and  rotation  of  the  body  (6  degrees  of  freedom). 


2.  Fundamentals  of  point  absorbers 

The  simulation  of  point  absorbers  is  a  wave-body  interaction 
topic.  It  requires  knowledge  of  body  dynamics  and  wave  theory. 
Before  we  review  the  modeling  methods  for  point  absorbers,  we 
review  the  fundamentals  of  body  dynamics  and  wave  theory. 


2.2.  Point  absorber  dynamics 

The  motion  of  a  point  absorber  in  six  degrees  of  freedom  is  plot¬ 
ted  in  Fig.  7.  The  dynamics  of  the  device  can  be  expressed  by  the 
following  equations, 

mbAt  =  Fhd  +  Fres  +  FpiO  +  Fext, 

„  ~  (1) 
^  x  IgS2  —  M frd  +  Mres  +  Mpro  +  Mext, 

where  mb  is  the  mass  of  the  body,  at  is  the  acceleration  vector 
for  the  translation,  £2  and  a^  are  the  angular  velocity  and  accel¬ 
eration  vectors,  lg  is  the  moment  of  inertia  tensor  at  the  center 
of  gravity,  F  and  M  are  the  resulting  force  and  moment  acting  on 
the  body.  The  subscripts  ‘/id’  and  ‘res’  represent  the  hydrodynamic 
excitation  component,  the  restoring  component,  respectively.  The 
subscripts  ‘PTO’  and  ‘ ext ’  represent  the  power  take  off  (PTO)  forces 
and  the  addition  constraints,  such  as  mooring  line  forces.  Faltin- 
sen  [35]  comprehensively  reviewed  the  wave-induced  motion  and 
loads  on  offshore  structures  due  to  linear,  nonlinear,  and  viscous 
effects.  Note  that  the  paper  is  mainly  concerned  with  the  prediction 
of  the  hydrodynamic  components  of  the  system.  However,  Evans 
[36]  and  Pizer  [37]  showed  that  additional  constraints  might  fur¬ 
ther  reduce  the  motion  response  of  the  device.  The  effect  of  the 
mooring  constraint  depends  on  its  configuration  and  design.  A  dis¬ 
cussion  about  the  influence  of  the  mooring  system  on  the  power 
generation  performance  of  a  two-body  wave  energy  system  was 
presented  in  Ref.  [38]. 

Generally,  single-body  point  absorber  converts  energy  from  the 
relative  motion  between  the  floating  buoy  and  the  seabed.  Two- 
body  system  extracts  energy  from  the  relative  motion  between  the 
two  oscillating  bodies,  a  float  section,  and  a  reaction  section.  The 
float  is  only  allowed  to  move  freely  in  one  degree  of  freedom,  with 
respect  to  the  reaction  section  (generally  in  heave).  When  the  PTO 
component  was  included,  the  additional  PTO  force  FPT0  was  con¬ 
sidered.  A  mass-spring-damper  system  can  be  used  to  represent 
the  PTO  mechanism,  and  PTO  force  becomes 


Fpro  =  -CpTOurel  ~  k zreh 


where  zre/  and  urei  are  the  relative  motion  and  velocity  of  the  two 
bodies,  k  is  the  spring  stiffness,  CPT0  is  the  power  absorption  (damp¬ 
ing)  coefficient.  The  generated  power  from  PTO  is  defined  as 

PpTO  =  CpTOUreP  (3) 

which  is  proportional  to  the  square  of  the  relative  translational 
velocity  of  the  two  sections.  To  optimize  the  power  generation, 
different  control  strategies  can  be  applied  to  the  system.  An  intro¬ 
duction  to  latching  phase-control  was  given  in  Ref.  [39],  and  a 
review  on  the  implementation  of  a  spring  in  phase  control,  which 
is  called  reactive  phase-control,  was  described  in  Ref.  [26]. 


2.2.  Free  surface  wave  theory 


The  free  surface  waves  are  a  representation  of  various  forces, 
e.g.,  the  wind  or  a  ship,  acting  on  and  deforming  the  fluid  surface 
against  the  action  of  gravity  and  surface  tension.  Naturally,  waves 
occur  in  all  sizes  and  forms.  Depending  on  their  magnitude,  they 
act  on  the  fluid  along  with  other  environmental  conditions,  such 
as  bottom  topography,  temperature  and  fluid  density.  In  general, 
waves  can  be  described  as  linear  or  nonlinear,  regular  or  irregular, 
unidirectional  or  omni-directional.  The  analytical  solutions  of  free 
surface  waves  were  derived  based  on  the  potential  flow  method, 
where  the  flow  is  assumed  to  be  incompressible  and  irrotational. 
The  velocity  potential  0(X,t)  therefore,  is  obtained  by  solving  the 
Laplace  equation 

V2</>  =  0,  (4) 


The  exact  solution  of  the  above  equation  is  very  complex.  We 
start  the  review  at  the  linear  wave  theory,  which  is  obtained  using 
simplified  assumptions.  Linear  wave  theory,  also  referred  to  as 
small  amplitude  wave  theory  and  airy  wave  theory  [40,41],  is  the 
simplest  solution  (first-order  approximation)  for  the  flow  field.  Lor 
progressive  gravity  waves  of  period  T,  amplitude  H,  and  water  depth 
d,  linear  wave  theory  assumes  FI  is  much  smaller  than  d  and  wave¬ 
length  X,  and  the  boundary  conditions  are  linearized.  The  kinematic 
free  surface  boundary  condition  (KLSBC)  and  the  dynamic  free  sur¬ 
face  boundary  condition  (DLSBC)  are  then  reduced  to 


KFSBC  :  ^  ^  =  0  at  Z  =  0, 

9z  dt 

DFSBC  :  -~$-+gV  =  0  at  z  =  0, 

ot 


where  g  is  gravity.  Both  the  linearized  KLSBC  and  DLSBC  are  defined 
at  the  mean  free  surface.  The  flow  properties  of  linear  waves  can 
be  expressed  as 

Wave  elevation:  rj  =  —  cos[/<(x  -  ct)],  (6) 


Velocity  potential  : 


<P  = 


gH  cosh [k{z  +  d)] 
2  co  cosh(kd) 


sin[k(x  -  ct)], 


Dispersion  relation  :  c2 


^tanh(kd), 


Horizontal  flow  velocity  : 


u  = 


ttH  cosh [k{z  +  d)] 
T  sinh  (kd) 


cos[/<(x  -  ct)], 


Vertical  flow  velocity : 


w  = 


ttH  sinh[/<(z  +  d)] 
T  sinh(kd) 


sin[/<(x  -  ct)], 


(10) 


where  x  and  z  are  the  horizontal  and  vertical  positions,  k  (=2 n/X) 
is  the  wave  number,  c  (=A/T)  is  the  wave  speed,  and  co  is 
the  wave  frequency.  Depending  on  the  water  depth,  the  waves 
can  be  categorized  as  shallow  water  waves  when  d/A<l/20,  as 
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Fig.  8.  The  applicability  of  various  wave  theories  [45,47]. 


intermediate-depth  waves  when  1  /20  <  d\ X  <  1  /2,  and  as  deep  water 
waves  when  d/ X  >  1  /2. 

When  the  wave  amplitude  is  not  small,  the  waves  become  non¬ 
linear.  The  linearized  KFSBC  and  DFSBC  are  not  satisfied,  and  the 
nonlinear  KFSBC  and  DFSBC  can  be  written  as 


KFSBC  : 


dr]  90  dr]  90 

dt  +  dx  dx  dz 


at  z  =  rj{x,  t), 


P  1 
DFSBC :  -  +  - 

P  2 


90 

9F+gz 


C(t) 


at  z  =  t]{x,  t ), 


(11) 


where  C(t)  is  constant. 

The  nonlinear  solution  can  only  be  obtained  by  solving  Eqs. 
(5)-(8).  The  nonlinear  higher  order  solution  is  generally  the  super¬ 
position  of  additional  components,  which  are  at  higher  wave 
frequencies,  on  the  fundamental  linear  wave-theory  terms.  There 
have  been  numerous  nonlinear  higher  order  wave  theories  devel¬ 
oped  since  the  1800s.  The  Stokes  theory  [42,43]  is  the  most  well 
known  and  studied.  An  extension  of  the  fifth-order  Stokes  waves 
was  presented  by  Skjelbreia  and  Hendrickson  [44],  and  it  has  been 
widely  used  in  engineering  applications.  The  wave  profile  and  the 
velocity  potential  of  the  fifth-order  Stokes  waves  are  given  as 

5 

Wave  elevation  :  kr]  =  y^An  cos[nk(x  -  ct)],  (12) 

n= 1 


Velocity  potential  :  —  =  cosh[n/<(z  +  d)]  sin [nk(x  -  ct)].(  13) 

n= 1 

The  applicability  of  various  wave  theories  was  studied  in 
Refs.  [45,46],  and  the  figure  plotted  by  Le  Mehaute  is  shown  in 
Fig.  8.  It  illustrates  the  validity  limitation  of  these  approximation 
approaches  for  different  wave  conditions.  A  review  of  the  introduc¬ 
tion  of  wave  theories  and  the  applications  to  wave  load  prediction 
on  offshore  structures  is  provided  by  Sarpkaya  [47],  and  a  more 
comprehensive  review  of  the  theories  and  the  applications  of  linear 
and  nonlinear  waves  is  given  in  Ref.  [48]. 

In  the  real  world,  waves  are  not  monochromatic.  There  are 
a  large  variety  of  waves  with  different  frequencies,  phases,  and 
amplitudes.  For  an  adequate  description  of  the  free  surface,  a  large 
number  of  waves  must  be  superimposed  together.  Usually,  we 
use  spectra  to  describe  them.  Furthermore,  the  behavior  of  waves 


is  rather  random.  In  this  situation,  the  waves  are  called  irregular 
waves.  The  random  wave  field  can  be  approximated  by  using  an 
infinite  sum  of  sinusoidal  components  propagating  with  various 
wavelengths.  Pierson-Moskowitz  spectrum  [49]  and  Jonswap 
spectrum  [50]  are  the  two  most  widely  used  spectra  profiles.  Refer 
to  textbooks  such  as  [48,51]  for  wave  theory  and  [35,47,52]  for 
wave-body  interactions. 


3.  Analytical  maximum  power  extraction  studies 


Analytical  methods  can  be  sorted,  with  or  without  detailed 
descriptions,  by  WEC  systems.  The  first  analytical  study  of  a  point 
absorber’s  power  output  focused  on  maximum  efficiency,  without 
the  detailed  description  of  the  device,  in  1975  [53,54].  The  lin¬ 
earized  equation  of  motion  for  the  point  absorber  device  in  heave 
is 


(m  +  mr)z  +  RyZ  +  Sz  —  Fpro  +  Text? 


(14) 


where  m  is  the  mass  of  the  body,  mr  is  the  hydrodynamic  coef¬ 
ficient  of  added  mass,  Rr  is  the  radiation  damping  coefficient,  S 
is  stiffness  of  hydrostatic  restoring  force,  FPT0  is  the  vertical  force 
due  to  power-take-off  mechanism,  and  Fext  is  excitation  force.  The 
maximum  energy  absorption  Pmax  for  different  incident  wave  fre¬ 
quencies  can  be  reached  when  the  oscillating  body  system  is  at 
resonance.  Based  on  the  linear  deep-water  wave  assumptions,  the 
time  averaged  optimal  power  absorption  at  resonance  is 

Pmax  —  7T5- I ^ext  P-  (15) 

O  l\j~ 

An  absorption  width  Lmax  is  defined  as 

T  Pmax  A 

Lmax=_  =  _  (16) 


pg2TH 2 
32 71 


(17) 


where  J  is  the  wave  energy  flux  for  linearized  deep  water  waves. 
Note  that  Evans  [55],  Mei  [56]  and  Newman  [57]  also  independently 
derived  the  same  result  in  Eq.  ( 1 6)  using  similar  approaches  in  1 976. 

The  maximum  power  that  can  be  absorbed  from  waves  and  the 
absorption  width  is  derived  based  on  linear  wave  theory  and  con¬ 
sideration  of  wave  radiation  and  diffraction.  Budal  and  Falnes  [58] 
proposed  another  upper  limit  (also  referred  to  as  Budal’s  upper 
bound)  for  wave  power  that  is  absorbed  by  a  given  submerged  body 
volume  V 


Pu  TtpgH 
V  <  4  T 


(18) 


where  the  diffraction  effect  is  neglected,  and  the  size  of  the  body 
is  assumed  to  be  small  compared  with  the  wave  length.  The  two 
curves  represent  the  theoretical  prediction  of  the  maximum  wave 
energy  that  can  be  captured  by  a  semi-submerged  heaving  body 
with  optimum  heaving  amplitude.  The  power  output  for  a  practical 
device  is  below  the  two  upper  limits,  and  the  maximum  power  that 
can  be  absorbed  is  about  half  of  Pmax,  when  a  phase  control  method 
is  applied  [29]. 


4.  Modeling  methods 

The  dynamic  response  of  point  absorbers  can  be  calculated 
by  solving  Eq.  (1),  and  power  generation  can  be  predicted  using 
Eq.  (3).  Theoretically,  the  hydrodynamic  forces  that  affect  the 
point  absorber  motion  can  be  obtained  analytically,1  empirically 


1  The  analytical  method  can  be  derived  by  an  equivalent  linearization  of  the 
empirical  calculation.  Therefore,  the  former  can  be  viewed  as  a  special  case  of  the 

latter. 
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or  numerically,  depending  on  the  modeling  purpose,  the  opera¬ 
tional  condition,  and  the  detailed  design.  These  modeling  methods 
evolved  from  modeling  ship  hydrodynamics  and  offshore  floating 
structures.  Therefore,  we  apply  the  methods  for  naval  architecture 
and  offshore  engineering  when  we  review  modeling  methods  for 
predicting  the  hydrodynamics  of  point  absorber  systems. 

4.2.  Analytical  methods 

A  common  approach  to  determine  the  hydrodynamic  force  ana¬ 
lytically  is  to  presume  that  the  force  is  the  sum  of  wave  diffraction 
and  wave  radiation  components, 

Fhd  =  Fr  +  Fd  ,  (19) 

where  the  subscripts  ‘R’  and  ‘D’  represent  the  radiation  and  diffrac¬ 
tion  terms,  respectively.  The  radiation  components  include  an 
added-mass  component  and  a  radiation  damping  term. 

Fr  =  AX  +  BX,  (20) 

where  A  and  B  are  the  hydrodynamic  added-mass  and  damping 
coefficients.  The  wave  excitation  force,  induced  by  wave  diffrac¬ 
tion,  includes  a  diffraction  component  Fdf,  due  the  present  of  the 
floating  body,  and  a  Froude-Krylov  force  component  Ffi<}  generated 
by  undisturbed  incident  waves. 

FD  =  f  ^nds=  j  +  ^nds  =  Ffk  +  Fdf,  (21) 

J Sw  J Sw 

where  Sw  is  the  immersed  surface  on  the  body,  n  is  the  unit  nor¬ 
mal  vector,  0/  is  the  incident  wave  velocity  potential  and  0D  is  the 
diffraction  velocity  potential. 

The  velocity  potential  are  then  solved  analytically,  and  most 
of  the  analytical  solutions  were  derived  through  the  use  of  linear 
wave  theory.  Many  analytical  studies  can  be  found  in  the  literature. 
For  examples,  Newman  [59]  calculated  the  excitation  force  and 
moment  though  the  use  of  Haskind  relation.  Yeung  [60]  proposed 
an  analytical  method  for  calculating  the  added  mass  and  damp¬ 
ing  coefficients  for  offshore  engineering  applications.  The  paper 
suggests  a  method  to  calculate  the  coefficient  using  the  Eigen  func¬ 
tions.  In  addition,  an  analytical  derivation  of  the  hydrodynamic 
coefficients  of  a  floating  hemisphere  undergoing  forced  periodic 
oscillations  was  presented  in  Ref.  [61  ].  The  analytical  methods  have 
also  been  widely  used  in  ocean  engineering  problems,  especially  for 
an  array  of  cylinders  [62,63]. 

In  reality,  the  viscous  damping  effect  is  important  for  modeling 
the  hydrodynamics  of  point  absorbers,  and  neglecting  the  viscous 
drag  will  result  in  an  unrealistic  large  prediction  on  the  response  of 
the  device  and  on  power  extraction  performance,  particularly  near 
the  resonance.  Discussion  on  energy  loss  due  to  viscous  damping 
can  be  found  in  [34,64].  A  common  approach  to  include  the  viscous 
drag  effect  in  the  modeling  is  to  add  a  quadratic  damping  term  to  the 
non-viscous  hydrodynamic  forces.  More  details  on  the  calculation 
of  viscous  damping  force  are  discussed  in  Section  4.3. 

In  the  past  decade,  many  researchers  have  used  the  analyt¬ 
ical  method  to  calculate  the  power  output  of  point  absorbers 
[26,65-67].  In  addition  to  a  single  unit,  this  method  has  been  used 
to  study  arrays  of  WEC  systems.  For  example,  through  the  use  of 
linear  theory  and  multiple-scale  approximation,  Garnaud  and  Mei 
[68]  presented  an  analytical  solution  to  the  wave  energy  extrac¬ 
tion  using  infinite  strips  of  buoys  and  a  circular  array.  The  study 
showed  that  an  array  of  heaving  wave-absorbing  buoys  could  have 
a  broader  range  of  wave  frequencies  than  a  single  large  buoy,  with 
good  energy-absorbing  efficiency.  In  general,  analytical  methods 
are  very  efficient  for  providing  a  quick  performance  estimation  of 
point  absorbers  with  simple  geometries. 


4.2.  Boundary  integral  equation  methods 

Boundary  integral  equation  method  (BIEM),  also  referred  to  as 
boundary  element  method,  is  an  advanced  potential  flow  method 
that  can  be  used  for  handling  more  complicated  geometries.  The 
governing  equation,  Eq.  (4),  is  formulated  in  a  boundary  integral 
equation  form,  which  is  derived  by  introducing  a  Green’s  function 
G(p,q),  where  p  and  q  represent  the  field  point  and  the  source  point, 
respectively.  The  boundary  integral  equation  obtained  through  the 
use  of  Green’s  third  identity  gives 

«(p)</>(p)  +  J  </>(p)C„(p,  q)  -  </>n(p)G(p,  q)  ds  =  0,  (22) 

where  S  indicates  the  boundaries  of  the  fluid  domain,  of(p)  is  the 
internal  angle  formed  at  the  boundaries,  Gn  =  vGn  and  0n  =  v0  n, 
with  n  being  the  normal  vector  at  q.  The  boundary  surfaces  are 
discretized  into  panels.  By  using  Dirichlet-  and  Neumann-type 
boundary  conditions,  the  potential  flow  field  can  be  obtained  by 
solving  the  resulting  system  of  linear  equations  numerically. 

After  the  potential  flow  field  is  obtained,  the  pressure  on  the 
body  surface  is  evaluated  using  Bernoulli’s  equation, 

p  =  -p(c/>t  +  iv(/>-V(/>  +  gz)  .  (23) 

The  forces  F  and  moment  M  acting  on  the  body  are  calculated 
by  integrating  the  pressure  on  the  immersed  body  surface, 

F=  /  pnds;  M=  /  pnxrds,  (24) 

JsB  JsB 

where  r  is  the  distance  vector  from  the  center  of  rotation.  The 
viscous  damping  coefficient  can  be  estimated  based  on  the  empir¬ 
ical  solution,  following  the  same  approach  used  in  the  analytical 
method. 

Regarding  the  simulation  of  the  interaction  between  waves  and 
floating  bodies,  two  types  of  BIEM  approaches  have  been  used 
for  modeling.  One  follows  a  weakly  nonlinear  theory  that  uses  a 
perturbation  expansion  for  the  solution  (the  theory  of  wave  radi¬ 
ation  and  diffraction)  and  specifies  the  boundary  conditions  at  the 
mean  free  surface  and  body-surface  through  the  use  of  Taylor  series 
expansion.  The  other  method  uses  a  fully  nonlinear  time-domain 
approach. 

The  weakly  nonlinear  approach  has  the  advantage  of  having  the 
same  coefficient  matrix  for  the  system  of  equations  solved  at  every 
time  step.  One  can  solve  the  problem  in  the  time  domain  [69]  or 
in  the  frequency  domain,  which  is  the  method  used  in  the  panel 
code  WAMIT  [70].  In  the  frequency  domain  calculation,  radiation, 
diffraction,  and  excitation  forces  are  calculated  through  the  use 
of  the  linear  superposition  principle,  which  is  limited  to  smaller 
amplitude  waves  (weak  nonlinear  waves)  [71].  The  complexity  of 
the  method  increases  significantly  for  higher  order  schemes.  Exam¬ 
ples  of  studies  that  used  this  weak  nonlinear  approach  for  wave 
radiation  and  diffraction  problems  can  be  found  in  Refs.  [72-78]. 
A  discussion  on  the  use  of  linear  and  second-order  approaches  for 
predicting  the  wave  excitation  of  floating  bodies  was  given  in  Ref. 
[79]. 

The  fully  nonlinear  time  domain  approach  for  wave  analy¬ 
sis  was  first  proposed  by  Longuet-Higgins  and  Cokelet,  using  a 
Mixed  Eulerian-Lagrangian  (MEL)-type  method  for  modeling  sur¬ 
face  waves  [80,81].  Fully  nonlinear  boundary  conditions  were 
applied  on  the  instantaneous  water  free  surface  and  body  surface, 
and  the  free  surface  was  computed  using  a  tracking  method.  With 
BIEM,  the  computational  domain  is  discretized  along  the  domain 
boundaries,  including  the  free  surface  and  the  body  surface.  There¬ 
fore,  the  coefficient  matrix  for  the  system  of  equations  should 
be  updated  at  every  time  step,  and  the  method  is  only  solved  in 
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time  domain.  A  high  order  Runge-Kutta  scheme  is  often  used  to 
update  the  time  varying  free  surface  and  body  surface.  The  non¬ 
linear  time-domain  approach  is  more  accurate  for  modeling  the 
highly  nonlinear  interaction  between  waves  and  floating  bodies. 
Using  this  approach,  researchers  successfully  predicted  the  nonlin¬ 
ear  irregular  wave  radiation  [82]  and  diffraction  [83]  of  vertical 
cylinders,  as  well  as  the  ventilating  entry  of  a  surface-piercing 
hydrofoil  [84]. 

Additionally,  modeling  floating  bodies  in  numerical  wave  tank 
requires  a  boundary  that  generates  waves  and  absorbs  reflected 
waves  at  the  same  time,  and  a  boundary  that  absorbs  outgo¬ 
ing  waves.  Over  the  years,  many  efforts  have  been  made  to 
develop  a  numerical  method  to  simulate  these  types  of  bound¬ 
aries,  without  creating  additional  numerical  disturbances.  A  review 
on  unbounded  domains,  with  various  types  of  artificial  boundary 
conditions  and  absorbing  methods,  was  given  in  [85].  Most  of  the 
studies  between  the  1970s  and  the  1990s  were  based  on  potential 
flow  theory.  In  particular,  a  sponge-layer  method  was  proposed  by 
Israeli  and  Orszag  [86],  in  which  an  artificial  damping  layer  (zone) 
was  implemented  to  attenuate  the  outgoing  wave  before  it  reached 
the  outer  boundary.  Similar  approaches  also  were  proposed  in  Refs. 
[87,88]. 

The  frequency  domain  BIEM  is  the  most  commonly  used  BIEM 
method  for  modeling  the  hydrodynamics  and  the  power  genera¬ 
tion  performance  of  point  absorber  systems.  The  method  is  used 
for  calculating  the  non-viscous  hydrodynamic  forces  described  in 
the  analytical  method.  McCabe  et  al.  [89]  demonstrated  an  exam¬ 
ple  of  modeling  the  responses  of  axisymmetric  bodies  in  irregular 
waves  using  WAMIT.  The  method  was  combined  with  a  Laplace 
transfer-function  formulation  to  model  the  dynamics  of  the  floating 
body  in  time  domain.  Given  that  the  BIEM  can  be  used  for  modeling 
more  complex  geometry,  the  method  has  been  widely  used  in  the 
studies  of  not  only  point  absorber  systems  but  also  a  wide  variety 
of  WEC  systems.  For  examples,  the  method  was  used  to  study  the 
power  generation  performance  of  floating  pitching  devices  [9,90],  a 
modified  version  of  the  Salter’s  Duck  by  Cruz  and  Salter  [91  ],  and  a 
bottom-hinged  surge-pitching  plate  [92].  Moreover,  a  recent  study 
utilized  the  frequency  domain  BIEM  to  perform  a  series  of  studies 
on  the  power  extraction  performance  of  eight  selected  WEC  sys¬ 
tems,  including  five  heaving  absorber  systems  with  different  design 
principles  [64].  The  study  also  analyzed  the  annual  power  gener¬ 
ation  performance  and  the  capture  efficiency  at  five  different  site 
locations.  Moreover,  several  researchers  have  used  the  method  in 
the  study  of  control  and  tuning  strategies  for  heaving  absorbers 
[93-96].  For  the  analysis  of  array  systems  for  point  absorbers,  the 
interaction  effect  could  be  considered  through  the  use  of  a  multiple¬ 
scattering  method  [67,97,98].  In  particular,  Vicente  et  al.  [98]  used 
the  method  to  analyze  the  dynamics  of  the  array  system  with  inter¬ 
body  connections  and  bottom  slack-mooring  lines. 


4.3.  Viscous  drag  calculation  for  potential  flow  solution 


Generally,  the  additional  viscous  damping  term  follows  the  drag 
term  in  the  Morison’s  equation.  The  Morison’s  equation  was  pro¬ 
posed  by  four  authors,  J.R.  Morison,  M.P.  O’Brien,  J.W.  Johnson  and 
S.A.  Schaaf,  at  the  University  of  California,  in  1 950  [99],  and  the  drag 
term  was  given  as 


1 

2 


pCdAc{u  -  ub)\u 


(25) 


where  Cd  is  the  quadratic  drag  (viscous  damping)  coefficient,  p  indi¬ 
cates  the  density  of  sea  water,  Ac  is  the  characteristic  areas,  u  is 
the  flow  velocity,  and  ub  is  the  velocity  of  the  body.  Generally,  the 
value  of  viscous  damping  coefficient  depends  on  the  body  geom¬ 
etry,  the  Reynolds  number  and  the  Keulegan-Carpenter  number. 


More  details  of  the  use  of  quadratic  drag  approach  and  the  appro¬ 
priate  values  for  the  empirical  coefficient  can  be  found  in  [47,100]. 

In  addition,  the  empirical  viscous  damping  coefficient  data  are 
limited  to  existing  simple  geometry.  For  practical  point  absorber 
geometry,  the  hydrodynamic  forces  may  have  to  be  evaluated  by 
conducting  additional  wave  tank  tests  [101]  or  prescribed  motion 
CFD  simulations  [102,103]  if  existing  data  is  not  available. 


4.4.  Navier-Stokes  equation  methods 


The  viscous  effects  of  boundary  layer  separation,  turbulence, 
wave  breaking,  and  overtopping  are  important  to  the  prediction  of 
hydrodynamic  forces  of  the  devices.  The  potential  flow  method  can¬ 
not  capture  them.  Therefore,  modelers  often  turn  to  a  fully  viscous 
solution  by  implementing  the  Navier-Stokes  equation  methods 
(NSEM).  The  velocity  vector  U  and  pressure  p  of  the  incompress¬ 
ible  flow  field  are  obtained  by  solving  the  continuity  equation  and 
the  momentum  equations,  which  are  given  as 


VU  =  0, 

p{ Ut  +  U  •  VU)  =  -  Vp  +  ¥b  +  V  •  T, 


(26) 


where  Ut  is  the  time  derivative  of  U,  ¥b  is  the  body  force  vector  (such 
as  gravity),  and  T  indicates  the  stress  tensor.  The  computational 
domain  can  be  discretized  using  a  finite  difference/volume/element 
method,  and  the  resulting  system  of  linear  equations  can  be  solved 
directly,  or  preferably,  through  iterations  using  an  appropriate  lin¬ 
ear  equation  system  solver.  A  well  known  Semi-Implicit  Method  for 
Pressure-Linked  Equations  (SIMPLE)  algorithm  [104]  and  its  vari¬ 
ants  have  been  widely  used  for  solving  Navier-Stokes  equations. 
More  details  on  the  numerical  discretization  methods,  as  well  as 
the  physical  phenomena,  are  given  in  Refs.  [105,106].  The  force 
and  moment  calculation  is  similar  to  that  applied  in  the  potential 
flow  theory,  except  that  the  frictional  force  also  is  included  in  the 
calculation. 

Basically,  using  the  Navier-Stokes  equation-based  method  for 
modeling  the  hydrodynamics  of  floating  bodies  in  waves  involves 
the  calculation  of  the  free  surface  in  a  numerical  wave  tank  and 
the  simulation  of  the  turbulent  flow.  In  general,  two  types  of 
approaches  have  been  applied  for  predicting  the  free  surface.  One 
is  the  surface  tracking  method,  which  treats  the  free  surface  as 
a  sharp  boundary  and  updates  it  with  time  [107].  However,  the 
method  does  not  have  the  ability  to  model  wave  breaking  and 
overtopping.  The  other  is  the  interface-capturing  method,  where 
the  simulation  is  performed  on  a  grid  that  includes  both  air  and 
water  phases.  Marker-and-Cell  (MAC)  method  [108],  Volume  Of 
Fluid  (VOF)  method  [109],  and  the  level  set  approach  [110,111] 
are  the  most  often  used  interface-capturing  methods.  In  particu¬ 
lar,  the  VOF  method  is  the  most  widely  used  for  commercial  and 
open  source  computational  fluid  dynamics  (CFD)  codes.  Review  of 
these  methods  and  their  application  to  wave  hydrodynamics  can 
be  found  in  [106,112].  The  artificial  damping  layer  method  also  is 
widely  used  in  the  Navier-Stokes  equation-based  approaches  for 
absorbing  the  outgoing  and  reflecting  waves  in  the  wave-body 
interaction  simulation.  For  example,  Lara  et  al.  [  1 1 3  ]  combined  the 
sponge-layer  method,  with  an  internal-wave  maker  method  [114], 
to  model  the  wave  generating-absorbing  boundary  and  simulated 
the  interaction  between  waves  and  permeable  structures. 

In  general,  four  types  of  numerical  methods  are  applied 
for  turbulent  flow  modeling,  including  direct  numerical  sim¬ 
ulations  (DNS)  [115],  large  eddy  simulations  (LES)  [116,117], 
Reynolds-averaged  Navier-Stokes  (RANS)  methods,  and  detached 
eddy  simulations  (DES).  The  corresponding  computational  cost  for 
DNS  and  LES  is  very  high  DNS,  and  the  two  models  are  often 
viewed  as  a  research  tool  in  the  study  of  turbulence,  particularly 
for  most  floating  body  dynamics  problems.  Therefore,  the  RANS 
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method  is  the  most  widely  used  scheme  for  modeling  turbulent 
flow.  The  unsteadiness  of  the  turbulence  is  time  averaged  and 
the  computation  of  turbulent  flow  relies  on  simplified  engineering 
approximations.  DES  is  a  hybrid  method  that  combines  the  RANS 
and  LES  methods.  It  handles  the  near  wall  region,  with  a  RANS- 
type  model,  and  treats  the  rest  of  the  flow  field  using  a  LES-type 
method.  Different  types  of  RANS  models  are  designed  for  different 
turbulent  flows.  The  physical  phenomena  and  the  basic  assump¬ 
tions  for  each  turbulence  model  were  comprehensively  reviewed 
by  Wilcox  [118]. 

Using  NSEM  for  modeling  floating  body  dynamics  often  involves 
the  handling  of  moving  boundary  and  mesh  deformation.  Sev¬ 
eral  methods  have  been  developed  over  the  years,  including: 
(1)  boundary-fitted  mesh  morphing  method,  which  is  the  most 
common  approach  and  typically  involves  the  use  of  Arbitrary 
Lagrangian-Eulerian  methods  [119];  (2)  overset  method  (also 
referred  to  as  Chimera  approach),  which  use  a  multi-mesh  system 
that  solve  the  problem  using  separately  generated  but  overlapping 
meshes  through  iterations  [120];  (3)  immersed  boundary  method, 
which  introduces  a  body  force  term  in  the  governing  equation  to 
represent  the  boundary  surface  [121];  (4)  Smooth  Particle  Hydro¬ 
dynamics  (SPH)  method  [122],  which  is  a  mesh  free  Lagrangian 
method  and  uses  a  smoothing  kernel  to  approximate  a  delta  func¬ 
tion,  with  many  sampling  points  in  the  fluid  domain.  Generally, 
boundary-fitted  mesh  morphing  method  and  the  overset  method 
are  the  most  widely  used  models  for  solving  wave  and  floating  body 
dynamics  problems. 

NSEM  has  often  been  used  for  studying  the  complex  nonlinear 
wave  and  floating  body  interaction  in  ocean  engineering  problems. 
For  example,  Moctar  et  al.  [123]  used  a  RANS  model  to  analyze  the 
wave  load  of  a  jack-up  platform  under  freak  waves,  where  the  VOF 
method  was  adopted  for  the  free-surface  calculation  and  was  cou¬ 
pled  with  a  finite  element  model  for  the  structure’s  stress  analysis. 
The  study  demonstrated  the  significance  of  using  a  Navier-Stokes 
type  approach,  where  the  effects  of  wave  run-up  on  the  platform 
leg  and  the  impact-related  wave  loads  on  the  hull  were  all  con¬ 
sidered.  Their  study  also  showed  that  the  difference  between  the 
shear  and  overturning  moment  for  the  platform  under  freak  waves 
calculated  from  the  Morison’s  approach,  and  those  obtained  from 
the  RANS  calculation,  can  increase  up  to  25%,  especially  for  larger 
and  breaking  waves. 

A  few  studies  on  the  fluid-structure  interaction  of  WEC  sys¬ 
tems  were  performed  using  the  Navier-Stokes-type  approaches. 
A  RANS  code  (COMET)  was  applied  by  Agamloh  et  al.  [124]  to 
model  the  dynamics  of  a  single  buoy  and  a  double  buoy  system 
in  waves.  The  PTO  mechanism  was  considered  as  the  damping 
system  based  on  the  energy  that  was  captured.  Studies  on  the 
Wave-driven,  Resonant,  Arcuate-action,  Surging  Point-Absorber 
(WRASPA)  were  presented  by  Bhinder  et  al.  [125,126].  The  authors 
performed  a  series  of  studies  using  a  RANS  model  to  investigate 


the  hydrodynamics  of  a  floating-point  absorber  [127].  The  power 
take  off  was  represented  using  a  mass-spring-damper  system,  and 
the  maximum  power  that  this  particular  device  generated,  with  a 
constant  value  spring-damper  system,  was  examined. 

5.  Discussion 

In  Section  4,  we  have  summarized  the  existing  numerical  meth¬ 
ods  for  modeling  the  dynamic  response  of  point  absorber  in  waves 
and  for  predicting  the  power  generation  of  the  system.  Here,  we 
provide  further  discussion  and  insights  on  these  methods. 

Table  1  illustrates  the  features  of  each  method  for  modeling 
point  absorber  systems.  Depending  on  the  purpose  of  the  study, 
each  method  has  its  advantages  and  limitations.  The  analytical  solu¬ 
tion  only  can  be  used  for  simple  geometries,  while  BIEM  and  NSEM 
can  be  used  for  more  complicated  and  realistic  ones.  The  analytical 
method  and  the  frequency  domain  BIEM  model  the  system  dynam¬ 
ics  in  a  sequential  approach  and  use  the  treatments  of  radiation 
and  diffraction  problems  for  estimating  the  hydrodynamic  forces. 
Therefore,  the  utilization  of  these  two  methods  is  generally  limited 
to  relatively  smaller  waves.  On  the  contrary,  the  time  domain  BIEM 
and  NSEM  fully  couples  the  flow  field  simulation  with  the  solution 
for  the  equation  of  motion  at  each  time  step.  The  time  domain  BIEM 
is  feasible  for  capturing  the  fully  nonlinear  interaction  between 
waves  and  the  floating  body,  except  for  wave  breaking,  overtopping 
and  viscous  flow  separation.  A  floating  body  dynamics  CFD  (NSEM) 
simulation  is  needed  to  fully  capture  all  the  nonlinear  interaction 
between  waves  and  the  device,  including  viscous  flow  separation, 
turbulence,  wave  breaking  and  overtopping. 

An  example  for  modeling  the  power  generation  performance  of 
a  two-body  floating-point  absorber  system  is  presented  to  demon¬ 
strate  the  utilization  of  these  methods.  The  two  most  commonly 
used  methods,  a  frequency-domain  BIEM  and  a  NSEM  (CFD),  were 
selected.  In  the  frequency  domain  BIEM,  the  hydrodynamic  forces 
were  calculated  using  WAMIT,  and  the  viscous  damping  coeffi¬ 
cient  was  estimated  based  on  the  empirical  solution  of  [47,100]. 
The  details  of  the  device  geometry  and  the  numerical  settings  in 
the  CFD  simulation  were  described  in  Ref.  [127].  To  compare  the 
results  obtained  from  the  two  methods,  the  computational  time  for 
each  model  is  shown  in  Table  2,  and  the  power  generation  perfor¬ 
mance  is  plotted  against  incident  wave  periods  in  Fig.  9.  It  is  found 
that  the  frequency-domain  BIEM  is  much  more  cost-effective.  How¬ 
ever,  the  selection  of  viscous  damping  coefficients  has  a  significant 
influence  on  the  frequency-domain  BIEM  results,  particularly  near 
resonance.  Note  that  the  maximum  power  generally  occurs  when 
point  absorber  systems  are  close  to  resonance.  The  calculation  of 
viscous  damping  is  important  to  predict  the  power  and  may  lead 
to  significant  uncertainty  in  the  prediction  results.  For  example, 
Babarit  et  al.  [64]  showed  that  if  the  drag  coefficients  were  varied 
by  a  factor  of  4,  the  effect  on  power  generation  could  be  up  to  30%. 


Table  1 

Features  of  hydrodynamic  modeling  methods. 


Methods 

Features 

Viscous  effects 

Dynamic  coupling 
complexity 

Wave/body  interaction 

Accuracy 

Hydrodynamics  modeling 

System  dynamics  modeling 

Analytical 

Empirical 

Radiation  and 

Linear 

Affected  by  the  prediction 

Limited  to  waves  with 

Frequency  domain 

Empirical 

diffraction  theory 
Radiation  and 

Non-breaking  waves 

of  viscous  damping 

Affected  by  the  prediction 

relatively  smaller 
amplitude 

BIEM 

Time  domain  BIEM 

Empirical 

diffraction  theory 

Fully  coupled  time 

Non-breaking  waves 

of  viscous  damping 

Affected  by  the  prediction 

Very  good 

NSEM:  floating-body 
dynamics  CFD 

Solves 

Navier-Stokes 

equation 

domain  simulation 

Fully  coupled  time 
domain  simulation 

Fully  nonlinear  and 
breaking  waves 

of  viscous  damping 

Very  good 

Fully  nonlinear 
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Table  2 


Run-time  benchmark  for  modeling  point  absorber  WEC  systems. 

Hydrodynamics  modeling 

Discretization 

Number  of  panels/meshes 

Wall  clock  time 

Frequency  domain  BIEM  (WAMIT)  &  quadratic  drag 
NSEM:  floating  body-dynamics  CFD 

Boundary  discretization 
Computational  domain 
discretization 

5096  (low-order  panels) 

million  (finite  volume  mesh) 

2804  sa  (30  wave  frequencies) 

12  h  on  64  coresb  (one  wave  frequency) 

a  On  a  3.33  GHz  Intel  i5  processor  with  8  GB  of  memory. 

b  Each  compute  node  consists  of  dual  socket/quad-core  2.93  GHz  Intel  Nehalem  processor,  with  12  GB  of  memory  shared  by  all  8  cores. 


For  operational  wave  conditions,  where  waves  are  relatively 
small  and  linear  theory  are  generally  applicable,  and  the  non- 
viscous  hydrodynamic  forces  are  dominated  by  radiation  and 
diffraction  components.  The  frequency  domain  BIEM  is  the  most 
commonly  used  method  for  modeling  the  power  generation  per¬ 
formance  of  point  absorber  systems.  However,  the  viscous  damping 
coefficients  need  to  be  carefully  selected.  The  drag  coefficient 
depends  on  the  device  geometry,  the  scale  and  the  relative 
velocity  between  the  body  and  the  flow  around  it.  The  drag  coef¬ 
ficient  becomes  much  larger  when  the  Reynolds  number  and 
the  Keulegan-Carpenter  number  are  smaller  [100].  In  addition, 
after  the  frequency  domain  hydrodynamic  forces  are  obtained,  the 
response  of  point  absorbers  can  be  calculated  using  a  frequency 
domain  approach  or  a  time  domain  approach.  The  frequency 
domain  solution  represents  the  steady-state  response  due  to  the 
sinusoidal  excitation  of  regular  waves.  The  time  domain  approach 
solves  the  Cummins  equations  [128]  and  are  often  used  for  mod¬ 
eling  the  dynamics  of  the  system  in  irregular  waves,  particularly 
for  the  study  of  optimal  control  strategies.  The  methods  and  the 
relationship  between  the  parameters  in  the  frequency  domain 
approach  and  those  in  the  time  domain  approach  were  well  sum¬ 
marized  by  Ogilvie  [129]. 

For  survivable  conditions,  a  time-domain  fully  nonlinear  BIEM 
model  or  NSEM  is  more  applicable.  In  reality,  point  absorbers  are 
designed  to  sustain  the  extreme  wave  load  during  storms.  The 
device  may  subject  to  waves  with  a  typical  100-year  significant 
wave  height  in  the  range  between  8  m  and  13  m.  Under  these  waves, 
breaking  and  overtopping  waves  often  occur  and  the  effect  of  the 
nonlinear  interaction  between  waves  and  the  floating  body  on  the 
dynamic  response  of  the  system  is  significant.  Only  NSEM  models 
have  the  capability  of  modeling  all  the  complex  nonlinear  hydrody¬ 
namics.  However,  longer  computational  time  is  required.  The  actual 
computational  time  generally  depends  on  many  factors,  including 


processor  speed,  convergence  criterion,  mesh  quality,  parallel  effi¬ 
ciency,  and  the  complexity  of  the  interaction  between  waves  and 
the  floating  body.  Thus,  one  may  have  to  compromise  accuracy 
against  computational  cost. 

6.  Conclusion 

Four  numerical  methods  for  modeling  the  point  absorber  sys¬ 
tems  have  been  reviewed  in  this  paper;  each  method  has  its  own 
strength  and  limitation.  For  idea  analysis,  analytical  method  is  suf¬ 
ficient  and  effective.  For  modeling  a  practical  point  absorber  WEC 
system,  frequency  domain  BIEM  and  NSEM  are  the  most  commonly 
used  method  depending  on  the  wave  condition.  In  operational 
condition,  the  frequency  domain  BIEM,  with  a  carefully  selected 
viscous  damping,  is  generally  used  for  power  generation  perfor¬ 
mance  analysis,  design  and  optimization.  They  tends  to  be  easier 
to  combine  with  the  modeling  techniques  for  other  aspects  of  the 
point  absorbers  such  as  the  power  system.  On  the  other  hand,  NSEM 
is  typically  used  for  detailed  analysis  or  for  survivability  analy¬ 
sis,  where  the  nonlinear  interaction  between  waves  and  the  WEC 
system  is  significant.  Time  domain  fully  nonlinear  BIEM  is  more 
appropriate  in  the  transition  between  these  those  conditions  and 
for  developing  active  control  strategies. 

In  additional  to  the  point  absorbers,  there  are  a  wide  variety  of 
WEC  systems  such  as  bottom  hinged  device  and  oscillating  water 
column.  Similar  modeling  philosophy  can  be  applied  to  them.  How¬ 
ever,  the  dominant  physical  phenomenon  for  each  system  varies. 
For  examples,  additional  nonlinear  effects  on  the  hydrodynamics 
of  the  system  due  to  wave  breaking  and  overtopping,  can  be  signif¬ 
icant  for  a  bottom-hinged  pitching  plate.  The  additional  constraint 
due  to  mooring  connections  can  be  essential  for  a  floating  device 
oscillating  water  column.  The  prediction  of  these  additional  forces 
and  the  level  of  complexity  of  the  interaction  between  forces  and 
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system  dynamics  will  affect  the  accuracy  of  the  modeling  method, 
particularly  the  frequency  domain  BIEM,  and  the  overall  computa¬ 
tional  time.  We  will  investigate  these  features  in  future. 
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